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Conventional methods of structural dynamic analysis are intuitive in 
nature. Equations of motion are obtained from an assumed model that, 
at best, has only a reasonable comparison to the actual structure. 
Because of their intuitive foundation, conventional analyses leave much 
to be desired with regard to the adequacy of proposed "fixes" or the 
dynamic effects of modifications such as gun pods, radar units, and 
external stores. 


This contract was initiated to develop a theory of structural dynamic 
testing which could be used to determine, directly from measurable test 
data, the equations of motion, eigenfunctions, and natural frequencies 
of a complex structure such as a helicopter. Within the framework of 
the idealized assumptions, the following major goals have been achieved: 


@ Theory derived and proven - an exact method for identifying the 
narameters in n equations of motion of an n-degr :e-of- freedom 
linear structure was developed and shown to be theoretically 
correct. 


@ Theory shown to be numerically sound - the method of implementing 
the theory was designed to eliminate ill-behaved matrices and 
excessive sensitivity to experimental or measurement error. 


@ Theory found to be experimentally practical - precisely controlled 
statistical computer experiments demonstrated that the theory is 
operable using measured input data of the type common to helicopter 
structural testing and with errors in excess of the accuracy of 
available testing equipment. 


A second contract has been awarded with the principal objectives of 
determining the adequacy of the chosen n degrees of freedom and the 
applicability of the theory of nonsymmetrical, three-dimensional structures. 
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ABSTRACT 


This volume contains a description of the computer program 
used to numerically test the system identification theory 

of Volume I. A sample computer run of an identification 
case which was discussed in Volume I is shown. The vari- 
ables used in the program are defined and operating instruc- 
tions for the program are presented in detail. 
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DESCRIPTION OF COMPUTER PROGRAM 


This complex computer program is designed to test and 
evaluate the methods derived in USAAVLABS Technical Report 
70-6, Volume I. The program is written in FORTRAN IV (E 
Level) for the IBM system 360. It has been run on a model 
40 using the Disk Operating System Version 3, Level l 
(release 17). It can be run without change on any IBM 360 
having 128,000 bytes of storage. It may be run on any com- 
puter having a FORTRAN IV compiler with only minor changes, 
providing sufficient storage is available. The program 
consists of a main program plus 13 FORTRAN subprograms. 


The program performs the following functions: (1) it simu- 
lates a sinusoidal frequency sweep test with errors from 
which are obtained the natural frequencies by iterating on 
the mobilities; (2) it simulates tests, as above, at the 
previously determined natural frequencies and carries out 
the identification of the parameters in the equations of 
motion; (3) it runs the above tests a number of times to 
determine the mean values and standard deviations in each 
of the identified parameters in order to evaluate the sen- 
sitivity of the identified parameters to measurement error; 
(4) after an identification has been performed, the frequency 
response of the identified system is calculated to assess 
the validity of the identified analytical model. 


For details of the input, logic, and output, see flow chart 
(Figures 1, 2, and 3),program listing, and sample computer 
runs. 


READ INPUT 
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Figure 2. Flow Chart of Computer Program. 


NO 


SET FREQUENCY TO 
NEXT NATURAL FREQ, 
YES 


SAMPLE STATISTICS 


OUTPUT 


REPETITIONS 
COMPLETED 


FORM DIAGONAL 


MASS MATRIX 


AVERAGE DAMPING CALCULATE MEANS 


FOR USE IN STANDARD DEVIATIONS 


POSSIBLE OF IDENT PARAMETERS 
SUCCEEDING 


RERUN 


OUTPUT 
GO TO NEXT CASE 


Figure 3. 


GO TO NEXT CASE 


Flow Chart. of Computer Program. 


OPERATING INSTRUCTIONS 


All input data must be in the following units. 


Mass - lb-sec*/in. 
Stiffness - lb/in. 
Frequencies - hz. 


The type of run is controlled by a constant (IC) in the 


first column of the heading card. 


has the following functions. 


IC DESCRIPTION OF RUN 


0 Frequency sweep 
Iteration for 2, yY 


1 End of run (last 
card of deck) 


2 Same as 0 
3 Same as 0 
4 Identification - 


statistical based 
on repeated tests 


INPUT 


Control card 
Mass, Stiffness 


Control card 
Mass, Influence 
Coef matrix 
from tape 


Control card 
only 

(Other data 

same as pre- 
vious case) 


Control card 
discrete 
forcing fre- 
quency 
(Other data 
same as pre- 
vious case) 


This control constant 


OUTPUT 


Exact natural 
frequencies 

and modes 
Driving point 
mobility 

Row of mobility 
matrix 

Iterated 2, yY 


Same as 0 


Same as 0 
except no exact 
frequencies 

and modes 


Mean and 
standard 
deviation 
parameters 
of equations 
of motion 


Ic DESCRIPTION OF RUN INPUT OUTPUT 


5 Identification - Same as 4 Iterated 2, y 
single test Accumulated 
parameters of 
equations of 
motion 
6 Rerun responses Control card Driving point 
following IC = 4 only mobility 
or 5 using iden- Row of mobility 
tified K, diagonal matrix 


M and average g 


Each case starts with a heading card, one case following 
another. The last card of the deck has a '1' in column l. 


The first case of any run must have IC = 0 or 2. 


DESCRIPTION OF INPUT CARDS 


Heading Card (every cage) 


Column 1 Ic Control constant (see above) 
2 IRI = 0 Iterates for 2, y during 
frequency sweep 
= 1 Omits iteration 
3 NAMP = 0 Response in complex velocity 
mobility 
= 1 Response in acceleration units, 
amplitude in g's and phase 
4-80 Head Arbitrary heading 


Control Card (every case) 


Column 1-2 ND No. of degrees of freedom 

6-10 G Scalar structural damping 
Not used for IC = 6 

11-15 NROW Row of mobility printed for IC = 0,2,3,6 
If NROW = 0, omit printing 
For IC = 4, NROW = no. of repetitions 

Not used for IC = 5 

16-20 NINC No. of sweep frequency increments (99 max) 
for IC = 0,2,3,6 
If NINC = 0, no sweep 
For IC = 4,5, NINC = no. of discrete 
frequencies 


eee | | 


Control Card 


Column 21-25 


26-30 


31-40 


41-50 


51-60 
61-70 


71-80 


Mass Cards (for IC 


(every case) - Continued 


OMO 


DOM 


OH 


PCT 


PCTB 


PHE 


IX 


M 


Used for IC = 0,2,3,6 only. OMO = first 
sweep frequency. If OMO = 0, .5 of the 
lowest natural frequency :.s used 

Used for IC = 0,2,3,6 only. DOM = sweep 
frequency increment. If DOM = 0, it is 
set equal to (1.25* highest natural 
frequency - OMO)/(NINC-1) 

High frequency used in iterating 

If OH = 0, then 1.25* highest natural 
frequency is used. 

Not used for IC = 6 

Random error on amplitude,error is 
uniformly distributed between -/+ PCT* 
each element of mobility. 

Not used for IC = 6 

Bias error on amplitude = PCTB* element. 
Not used for IC = 6 

Random phase error in degrees, distributed 
uniformly between -/+ PHE degrees. 

Not used for IC = 6 

Seed for random number generator con- 
trolling measurement errors. 

Not used for IC = 6 


= 0,2) 
Diagonal masses only. 8 per card 


(8F10.0) 
Use as many cards as necessary 


Stiffness Cards (for IC = 0) 


K 


Stiffness matrix. Lower triangle input 
only. Start each row on new card and 
end with diagonal element. 8 per card 
(8E10.0) 

Matrix will be symmetrized after input 


C Matrix from tape (for IC = 2) 


Unformatted record contains heading (20 words, 
first character blank); NX (order of matrix); 
((C(r,3), I = 1,NX), J = 1,NX) 


ee eee a te 


Discrete Frequency Caras (for IC = 4,5) 


Should contain NINC frequencies, each of 
which should converge on a4 different mode 
during iteration. 

8 per card (8F10.9) 


B(20,21) 
BF (20,21) 
BS (20,21) 


C (20,21) 
DUM (20) 
DPI (100,20) 


DPR(100,20) 
FRE (20) 


FREK (20) 


Gl 
GKS (20) 
GM (20) 


GMK (20) 
GMS (20) 


IT (100) 
ITN 
K (20,21) 


L 

MF (20) 

OM 

PHI (20,21) 
PHIK (20,21) 
SD (20,21) 
SG (20,21) 
SGF (20, 21) 
SGS (20,21) 


SK (20,21) 


SKS (20, 21) 
SM (20,21) 
SMS (20,21) 
STAR (10) 

TI (100,20) 


TR (100,20) 


LIST OF OUTPUT VARIABLES 


Identified influence coefficient matr 2x 

Mean identified influence coefficient matrix 
Standard deviation of identified influence 
coefficient matrix 

Exact influence coefficient matrix 

Iterated gamma (also other temporary uses) 
Driving point mobility (imaginary or phase 
angle) 

Driving point mobility (real or amplitude 

in g's) 

Calculated natural frequencies from inf. uence 
coefficient matrix 

Calculated natural frequencies from stiffness 
matrix 

Average identified damping coefficient 
Standard deviation of sums of columns of 
stiffness matrix 

Exact generalized mass from influence 
coefficient matrix 

Exact generalized mass from stiffness matrix 
Standard deviation of sums of columns of 
mass matrix 

Number of iterations 

Number of iterations 

Exact stiffness matrix 

Dummy index, also number of increment 

Mean the sums of the columns of mass matrix 
Forcing frequency 

Mode shapes from influence coefficient matrix 
Mode shapes from stiffness matrix 

Identified damping matrix 

Identified damping coefficient matrix 

Mean damping coefficient matrix 

Standard deviation of damping coefficient 
matrix 

Mean stiffness matrix or identified 
matrix 

Standard deviation of stiffness matrix 

Mean mass matrix or identified mass matrix 
Standard deviation of mass matrix 
Generalized parameter from gamma iteration 
Row of mobility matrix from sweep (imaginary 
or phase) 

Row of mobility matrix from sweep (real or 
amplitude in g's) 


stiffness 


A(20,21) 
AREP 

AREP1 

ARG 

CK 

CON 

GAM (100,20) 
I 


Il 


Icl 
IC2 


ND1 


PMAX 
$1,S2...S6 


VAL 
YI (20,21) 


YIH (20,21) 
YR(20,21) 
YRH (20,21) 
21 (20,21) 
ZR(20,21) 


LIST OF INTERNAL VARIABLES 


Dummy matrix used for cM, M /K 


= NREP 

= AREP-1 

Dummy variable 

Generalized stiffness 

Dummy 

Iterated gamma vectors (output through DUM) 
Dummy index 

Dummy index 

zIc+i1 

No. of input matrices read from tape 
(controls initial rewind) 

Dummy index 

Count of repetitions (1 to NREP) 

Dummy index 

Mean sum of columns of stiffness matrix 
Dummy seed for random numbers 

Square form of mass matrix 

Dummy variable 

Dummy variable 

Dummy variable 

Dummy variable 

No. of repetitions 

Forcing frequency squared 

Not used 

Last forcing frequency 

Maximum amplitude of modal elements 

Same as STAR(1), STAR(16) for second 
dominant mode 

Eigenvalue of gamma iteration (not used) 
Simulated measured mobility matrix (imaginary) 
at forcing frequency 

Simulated measured mobility matrix (imaginary) 
at high frequency 

Simulated measured mobility matrix (real) 
at forcing frequency 

Simulated measured mobility matrix (real) 
at high frequency 

Exact impedance (imaginary) at forcing 
frequency 

Exact impedance (real) at forcing frequency 
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LIST OF FORTRAN SUBROUTINES 


Converts mobility from velocity units to 
acceleration aS amplitude (in g's) and 
phase angle (in degrees) 

Complex inverse of complex matrix 


Incorporates measurement errors into 
simulated measurements 


Generalized function of form flag where 
f is a vector and A iS a Square matrix 


Inverse of a matrix 


Matrix iteration for eigenvalues and 
eigenvectors 


More general iteration on product of two 
matrices. Used for gamma iteration 


Matrix multiplication 

Caiculates complex impedance and mobility 
Special output for square matrix 

Random number generator 

Removes rows and columns from matrix 


Special matrix output 


tt 


COMPILATION AND SAMPLE RUN 


The "measured" mobility matrices [Y,] are generated, over a 
range of frequencies w, from exact [m], [k], and [d]. They 
are then polluted with known error, and values of [m], [k] 
and [dad] are identified, using the theory described herein, 
from the polluted mobilities. The cycle is completed by 
generating a third set of mobilities from the identified 
parameters. 


The sample run below refers to 8 percent random error and 

8 percent bias error in the mobility amplitude and 2° random 
error in the mobility phase. All three phases, generation of 
polluted [Y,] identification of parameters and rerun of 
identified model to regenerate [Y,], are included and are 
Signified, respectively, by IC = 2,5,6. The complete driving 
point response Y;;, and the third row, corresponding to the 
pilot seat of the transfer point response, Yij, are both 
presented over a frequency range from 1 Hz to 1400 Hz. 
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